Accepted for publication in the Astronomical Journal 



Mixed-Morphology Supernova Remnants in X-rays: Isothermal Plasma in 
HB21 and Probable Oxygen-Rich Ejecta in CTB 1 

Thomas G. Pannuti 1 

Spitzer Science Center, California Institute of Technology, MS 220-6, Pasadena, CA 91125; 
Current Address: Space Science Center, Department of Earth and Space Sciences, 235 Martindale 
Drive, Room 212E, Morehead State University, Morehead, KY 40351 

t . pannutiOmoreheadstate . edu 

Jeonghee Rho 

Spitzer Science Center, California Institute of Technology, MS 220-6, Pasadena, CA 91125; 
Current Address: SOFIA Science and Mission Operations /USR A, NASA Ames Research Center, 

MS 211-3, Moffett Field, CA 94035 

jrho@sofia.usra. edu 

Kazimierz J. Borkowski 
Department of Physics, North Carolina State University, Box 8202, Raleigh, NC, 27695 

kborkow@unity . ncsu . edu 
and 

P. Brian Cameron 

Division of Physics, Mathematics and Astronomy, 105-24, California Institute of Technology, 

Pasadena, CA 91125 

pbcOastro . caltech . edu 
ABSTRACT 

We present an analysis of X-ray observations of the Galactic supernova remnants 
(SNRs) HB21 (G89.0+4.7) and CTB 1 (G116. 9+0.2), two well-known members of the 
class of mixed-morphology (MM) SNRs. Our analysis draws upon observations of both 
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SNRs made with the Advanced Satellite for Cosmology and Astrophysics (ASCA): we 
have also used an archived Chandra observation of CTB 1 as part of this effort. We 
find a marked contrast between the X-ray properties of HB21 and CTB 1: in the case 
of HB21, the extracted spectra of the northwest and southeast regions of the X-ray 
emitting plasma associated with the SNR can be fit with a single thermal model with 
marginally enhanced silicon and sulfur abundances. For both of these regions, the 
derived column density and temperature are Nh ~ 0.3xl0 22 cm -2 and kT ~ 0.7 keV, 
respectively. No significant spatial differences in temperature or elemental abundances 
between the two regions are detected and the X-ray-emitting plasma for both regions is 
close to ionization equilibrium. Our Chandra spectral analysis of CTB 1 reveals that this 
source is likely an oxygen-rich SNR with enhanced abundances of oxygen and neon: this 
result is quite surprising for an evolved SNR like CTB 1. The high angular resolution 
Chandra observation of CTB 1 reveals spectral variations across this SNR: in particular, 
we have detected localized hard emission with an angular extent of ~ 1'. The extracted 
ASCA spectra for both the southwest and northeastern regions of CTB 1 cannot be fit 
with a single thermal component and instead an additional component is required to 
account for the presence of excess emission seen at higher energies. Based on our fits 
to the extracted ASCA spectra, we derive a column density Njj ~ 0.6 x 10 22 cm" 2 and 
a temperature for the soft thermal component of kT so j t ~ 0.28 keV for both regions. 
The hard emission from the southwest region may be modeled with either a thermal 
component with a temperature kT^ard ~ 3 keV or by a power law component with a 
photon index T ~ 2-3; for the northeast region the hard emission may be modeled with a 
power law component with a photon index T = 1.4. The detection of center-filled ejecta- 
dominated X-ray emission from HB21 and CTB 1 as well as other MM SNRs suggests 
a new scenario for the origin of contrasting X-ray and radio morphologies of this class 
of sources. Lastly, we have analyzed the properties of the discrete hard X-ray source 
1WGA J0001. 4+6229 which is seen in projection just inside the northeastern shell of 
CTB 1. Our extracted ASCA GIS spectra of this source are best fit using a power-law 
model with a photon index r=2.2^ t j ) '2 : this slope is typical for featureless power-law 
continua produced by rotation-powered pulsars. This source may be a neutron star 
associated with CTB 1. We find marginal evidence for X-ray pulsations from this 
source with a period of 47.6154 milliseconds. A deep radio observation of this source 
failed to reveal any pulsations. 

Subject headings: supernova remnants: general - supernova remnants: individual (HB21 
(G89.0+4.7), CTB 1 (G116.9+0.2), 1WGA J0001.4+6229) - X-rays: ISM 
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Introduction 



A new morphological class of supernova remnants (SNRs) known as the mixed-morphology 



1998 



Shelton et al 



(MM ) SNRs has been firmly established in the recent literature (|Rho fe Petre 
19991 ) . The defining characteristics of SNRs of this class include a shell- like radio morphology 



combined with a centrally-filled X-ray morphology. X-ray observations of these SNRs made with 
the Roentgensatellit (ROSAT), the Advanced Satellite for Cosmology and Astrophysics (ASCA), 
Chandra and XMM-Newton have found that the central X-ray emission from these SNRs is not 
non-thermal emission (as would be expected from a central plerion) but instead thermal from 
shock-heated swept-up interstellar material. Examp les of well-known MM SNRs include W28 
dRho fc Borkowskilbood). G290.1-0.8 (MSH 11-61 A) (jsiane et alibood ) and IC 443 (|Kawasaki et al 
20021 ). iRho &: Petrel (|1998l ) suggested that as many as 25% of the entire population of Galactic SNRs 



may belong to th is morphological class. Based on CO and infrared observations (|Koo et al .1 12001 



Reach et al.ll2005l ). it appears that many MM SNRs are interacting with nearby molecular and HI 
clouds. This result suggests a connection between the contrasting X-ray and radio morphologies of 
these SNRs with the interaction between these sources and adjacent clouds, but such a connection 
lacks a detailed theoretical basis at this time. Two leading scenarios have been advanced to explain 
the origin of the contrasting radio and X-ray morphologies of MM SNRs: in the first scenario - 
known as the evaporating clouds scenario - molecular clouds overrun by the expanding shock front 
of the SNR survive passage through the shock and eventually evaporate, providin g a source of ma- 
teria l that increases the de nsity of the interior X-ray-emitting plasma of the SNR (jCowie fc McKee 

In the second scenario 



1977 



White fc Long 



1991 



known as the radiative shell model - the SNR 
has evolved to an advanced stage where the shock temperature is low and very soft X-ray emission 
from the shell is absorbed by the interst ellar medium (ISM): therefore, the only d etectable X-ray 



emission is from the interior of the SNR (Cox et al. 



1999 



Shelton et al 



1999 



paper we analyze and discuss X-ray emission from tw o Galactic SN Rs 



20041 ). In the current 
HB21 and CTB 1 - which 



have both been previously classified as MM SNRs by lRho &; Petrel (| 19981 ). 



HB21 (G89.0+4.7) was discovered in a radio surve y by iBrown &: Hazard! (|l953l ). The radio 
angular extent of this SNR is large - 120x90 arcminutes (|Greenll2009al ) - and the radio morphology 
is a closed shell. The shell appears to be flattened along the eastern boundary and features bright 
regions along the northern and southern boundaries with a prominent indentation seen along the 
northern boundary. The radio emission from HB 21 is strongly polarized (3.7%=fc0.4%) with a pro - 
jected magnetic field tangential to the shell (|Kundulll97ll : iKundu et al.lll973l : iKothes et al.ll2006l ). 
suggesting that the shell was compressed during the radiative evol utionary stage of the SNR. The 
measured radio spectral index for this SNR is a ~ 0.4 (S u oc v~ a ) (ILeahyi 120061 : lGreenll2009bl ) but 
significant variations in the values (from 0.0-0.8 with a standard d eviation of 0.16) of the spectral 
index across the f ace of the SNR were observed by ILeahyi (|2006l ). Based on IRAS observations, 
Saken et al.l ([1992]) detected clumpy infrared filaments associated with this SNR. Filamentary op- 
tical emission from this SN R with an angular extent comparable to the radio shell was detected by 
Mavromatakis et al.l (|2007l ). Extensive evidence exists that indicates HB21 is interacting with ad- 
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jacent molecular cloud s: this evidence includes CO observations (jTatematsu et alJll99d : iKoo et al 
200ll ; lBvun et al.ll2006l ) as well as near- and mid- infrared observations o f the sites of shock- mole cular 
cloud interactions along th e northern and southern parts of the SNR (jShinn et al.ll2009i . |2010| ). HI 
observations toward HB21 (jTatematsu et al.lll990l ; lKoo fe Heileslll99lh have revealed a high velocity 
expanding shell associated with this SNR. Fi nally, HB21 ha s been the s ubject of prior pointed X-ra y 
observations made by Einstein I Leahy 1987 ) and ROSAT ( Rho 1995 1: Leahy Sz Aschenbach 1996 ): 
we include in this paper the ROSAT images presented previously by Rho ( jl995l ) in her PhD thesis 
work. No pulsars or 7-ray sources are believed to be associated with H B21: radio searches for a 
pulsar associated with this SNR were conducted by Biggs <fe Lvne ( 1996 ) and Lorimer et al. ( 19981^) 



but no candidate sources were found. The distance to HB21 is not well known: 



Tatematsu et al 



(|l99d ) argued for a distance of only 0.8 kpc based on an association betwe en the SNR and molecular 



mater ial that belongs to the Cygnus OB7 association ([Humphreys! Il978l ). However. lYoshita et al 



(|200ll ) suggested a distance of >1.6 kpc based on a correlation t hat t hose authors found between 
X-ray absorbing column density and extinction, and iBvun et all (|2006l ) suggested a distance of 1.7 
kpc based on CO observations. In this paper, we have adopted a distance of 1.7 kpc to HB21. 

CTB 1 (G116.9+0.2) was discovered in a survey of Galactic radio e mission at 960 MHz 
by IWilson &: Boltonl (119601) . Subsequent radio observations of this SNR (IVelusamv fc Kundu 



1974: 



Kothes et al 



Angcrho fer et al.lll977l : lLandecker et al 



1982 



Yar-Uvaniker et al. 



2004 



Tian &: Leahy 



2006 



20061 ) reveal a radio morphology that m ay be describe d as a nearly complete circular 



shell with a diameter of approximately 34 arcminutes (|Greenll2009bl ). The radio emission is bright- 
est along the western rim and a prominent gap is seen along the northern and northeastern sector 
of the circular emission. Like HB21, the magnetic field is aligned in the tangential direction, also 
suggesting that the shell was compressed during the radiative stage of the evolutio n of the SNR 



sec 



Kothes et al 



but co mpared to HB21 the degree of polarization is much lower (0.4%±0.1% 
(12006)). The measured spectral index of the observed rad io emission is a ~ 0.6 (lLandecker et al 



19821 : iKothes et alJbood : iTian fc Leahvlbood : lGreenll2009bh . CTB 1 has also been detected at opti 
cal wavelengths, in emission lines such as [ OIII] A5007 and [ SII] AA 6716, 6731; the observed optical 
shell-like m o rphol ogy closely matches the radio shell. The optical images of CTB 1 presented by 
Fesen et al.l (|1997l ) depict a remarkable contrast in the emission-line properties of this SNR: while 
[SII] emission is seen from roughly the entire optical shell (with the greatest amount of emission 
in the south), the [ OIII] emission appears to be almost entirely localized to the western rim of the 
shell. ISaken et al.l (119921 ) detected infrared emission from CTB 1 in the 60 /jm and 100 fim IRAS 
bands: an arc of infrared emission was seen in the 60 /im/100 /im ratio ima ge that appears to be 



coincident with the radio sh ell. CTB 1 was observed in X-rays with ROSAT (jHailey &: Craig 



Rho 


1995; 


Craig et al. 


1997) 



1995 



Craig et al.lll997l ): like other MM SNRs, the X-ray emission from this SNR (which has a 
thermal origin) lies interior to the radio and optical shells. Remarkably, the X-ray emission is also 
seen to extend through the known northern gap of the SNR. Like HB21, no pulsars or 7-ray s ources 
are believed to be associated with CTB 1: a radio search conducted bv lLorimer et al.l ()1998l ) for a 
pulsar revealed no candidate sources. Published distance estimates for this SNR have ranged from 
1.6 to 3.5 kpc; in this paper we have adopted a distance to CTB 1 of 3.1±0.4 kpc as measured by 
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Hailev fc Craid (| 1994 1. 



The organization of this paper is as follows: in Section [2] we describe the ASCA X-ray ob- 
servations of HB21 and CTB 1 and the Chandra observations of CTB 1, including details of data 
reduction. ROSAT and radio observations of these SNRs are also described in this section. In 
Section 3 we present the results of our spectral analyses for both HB21 and CTB 1 (in Section 
3.1 and Section 3.2, respectively). In Section 4 we discuss the nature of the hard discrete X-ray 
source 1WGA J0001. 4+6229: we have discovered weak evidence for pulsed X-ray emission from 
this source (which is seen in projection against CTB 1) and consider the possibility that it is a 
neutron star associated with CTB 1. We also present a search for radio pulsations from 1WGA 
J0001. 4+6229. Interpretations of our X-ray results for HB21 and CTB 1 are presented in Section 
5 and Section 6, respectively. We also detected hard X-ray emission from CTB 1: we discuss the 
nature of this emissio n in Section Our preliminary results of this paper have been presented in 



Pannuti Sz Rhol (12004) a fter which we note that similiar data sets were analyzed and presented by 



Lazendic fc Sland (120061). O ur primary results of HB21 are in agreement with and strengthen those 



of lLazendic fc Sland ()2006l ): for CTB 1, our paper presents extensive and thorough analysis of the 
Chandra and ASCA data in smaller-scale regions. We also report important new results for this 
SNR, including the probable detection of oxygen-rich ejecta from CTB 1 as well as spectral varia- 
tions across the object. In addition, an X-ray hard point-like source is identified and an analysis of 
its X-ray and radio properties is presented. Finally, the conclusions of this work are summarized 
in Section 8. 



2. Observations and Data Reduction 
2.1. ASCA Observations of HB21 and CTB 1 

Because the X-ray emission from both HB21 and CTB 1 cov er a large angular e xtent on the 



sky, two pointed observations were made of each SNR with ASCA (ITanaka et al.lll994l ). namely the 
southeast and northwest regions of HB21 and the southwest and northeast regions of CTB 1 (see 
Table[T] for details of these observations). These observations provided almost a complete spatial 
coverage of the X-ray emitting gas in both SNRs. The data reduction was conducted using the 
"XSELECT" program (Version 2.2), which is available from the High Energy Astrophysics Science 
Archive Research Center (HEASARC0). There were two types of instruments onboard ASCA - the 
Gas Imaging Spectrometer (GIS) and the Solid-State Imaging Spectrometer (SIS) - and both of 
these instruments were composed of two units denoted as GIS2, GIS3, SISO and SIS1, respectively. 
A single GIS unit sampled a field of view ~ 50' in diameter and a background count rate of the 
GIS is 5 x 10~ 4 counts cm -2 sec -1 keV -1 ; in comparison, a single SIS unit sampled a field of view 
approximately 44' x 44' in size. The nominal FWHM angular resolution of both the GIS and SIS 



http: / /heasarc. gsfc.nasa.gov 
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units were approximately 1 arcminute. The standard REV2 screening criteria were applied when 
reducing both the raw GIS and SIS datasets. We used "FMOSAIC" from the FTOOLS software 
package to generate an X-ray map by combining the GIS2 and GIS3 maps. We used the FTOOL 
"MKGISBGD" to prepare blank-sky background spectra and images for each extracted GIS source 
spectra: for these background datasets, point sources which are brighter than approximately 10~ 13 
ergs cm -2 sec -1 have been removed. Similarly, background spectra were generated using standard 
SIS blank-sky datasets for analyzing the extracted SIS source spectra. The standard GIS2 and 
GIS3 response matrix files (RMFs) were used for analyzing extracted the GIS source spectra while 
the FTOOL "sisrmg" was used to prepare RMFs for the extracted SIS source spectra. Finally, the 
FTOOL "ASCAARF" was used to prepare ancillary response files (ARFs) for each extracted GIS 
and SIS source spectra. 



2.2. Chandra Observation of CTB 1 

We have also a nalyzed an archival da taset from an additional X-ray observation of CTB 1 



made with Chandra (jWeisskopf et al.l l2002l ) . The corresponding ObsID of this observation is 2810 
(PI: S. Kulkarni) and it was conducted as part of a search for central X-ray sources associated with 
Galactic SNRs. This observation was conducted in FAINT Mode on 14 September 2002 with the 
Advanced CCD Imaging Spectrometer (ACIS) at a focal plane temperature of — 120°C such that the 
ACIS-I array of chips sampled a significant portion of the X-ray emitting plasma located interior to 
the radio shell of the SNR. The ACIS-I array is composed of four front-illuminated CCD chips: each 
chip is 8'.3 x 8'.3 and the field of view of the entire array is approximately 17' x 17'. These chips 
are nominally sensitive to photons in the 0.2 through 10 keV energy range: the maximum effective 
collecting area for each chip is approximately 525 cm 2 at 1.5 keV. The full-width at half-maximum 
(FWHM) angular resolution of each chip at 1 keV is 1" and finally the spectral resolution at 1 keV 
of each chip is 56. These data were reduced using the Chandra Interactive Analysis of Observations 
(CIAcI) package (Version 3.1) with the calibration database (CALDB) version 2.29. Standard 
processing was applied to this dataset: in particular, the task "acis_process_events" was used to 
generate a new event file where corrections for charge transfer inefficiency and time-dependent gain 
were applied. The data were also filtered for bad pixels, background flare activity and events which 
had a GRADE value of 1, 5 or 7. Finally, we applied the good time interval (GTI) file supplied by 
the pipeline (as well as the GTI file prepared when filtering for background flares) and the resulting 
total effective exposure time of the observation was 48.9 kiloseconds. 



Discrete sources were identified with the CI AO wavelet detection routine "wavdetect" (Freeman et al 



20021 ): in making a final image, these sources were excluded and the image was exposure-corrected 
and smoothed with the CI AO task "csmooth." We extracted spectra from several regions of the 
diffuse X-ray emission seen in the Chandra images using the CI AO task "dmextract": results of the 



2 http://cxc. harvard.edu/ciao/ 
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spectral analysis are presented in Section T3.21 When extracting spectra, we excluded point sources 
identified by "wavdetect" to help reduce confusion with emission from background sources. We 
prepared ARFs and RMFs using the CIAO tools "mkwarf ' and "mkrmf," respectively; background 
spectra were generated using a reprojected blank sky observation made with the ACIS-I array and 
available from the Chandra X-ray Center via the World Wide Web@ 



2.3. Additional Observations 



We also included ROSAT Position Sensitive Proportional Counter (PSPC) observations of 
H B21 i n our analysis: these observations were discussed already in some detail in the Ph.D thesis 
of iRhol (|1995l ). They extend beyond the two fields observed by ASCA and provide the complete 
spatial coverage of the SNR. The PSPC images are exposure and particle backgroun d corrected and 



merge d together using the analysis techniques for extended objects developed by ISnowden et al 



(|1994l ). The smoothing technique includes neighboring pixels within a circle of increasing radius 
until a selected number of counts is reached to optimize the signal-to-noise ratio. 

Lastly, we have augmented the X-ray datasets analyzed in this w ork with radio data provided 
by the Canadian Galactic Plane Survey (CGPS) (jTavlor et al.ll2003l ). From this survey we have 
obtained radio images of HB21 and CTB 1 at the frequencies of 408 MHz and 1420 MHz. The 
angular resolution and sensitivity of the 408 MHz data are 3/4 x 3/4 esc 5 and 0.75 sin 5 K (3.0 mJy 
beam" 1 ), respectively, while the angular resolution and sensitivity of the 1420 MHz data are 1' x 1 ' 
esc 5 and 71 sin 5 K (0.3 mJy beam -1 ), respectively. The reader is referred to lTavlor et al .1 ( 20031 ) 
for more description about the CGPS radio observations and the accompanying data reduction 
process. 



3. Results 
3.1. HB21 

In Figure[l]we present our broadband (0.7-10.0 keV) exposure-corrected mosaicked ASCA GIS 
image of HB21: we have overlaid radio emission contours using CGPS observations at a frequency 
of 408 MHz to show the extent of the SNR radio shell. In Figure [2] we present a mosaicked 
ROSAT PSPC image of HB21 (with the same radio contours overlaid) which depicts the entire 
extent of X-ray emission from the SNR (which extends beyond the two fields observed by ASCA). 
It is clear from inspection of these images that the X-ray emission is located in the interior of 
the well-defined SNR radio shell: this combination of X-ray and radio morphologies exemplifies the 
defining characteristics of mixed-morphology SNRs. The bulk of the interior X-ray emitting plasma 



3 See |http://cxc.harvard.edu/contrib/maxim/acisbg/| 
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is located j ust s outh of a prominent bend in the northern edge of the radio shell. We note that 



Koo et al.l (|200ll ) detected broad CO emission lines from the location of this bend and lShinn et al 



(|2009i ) presented near- and mid-infrared images of this same region which showed shock-cloud 
interaction features. Both of these studies indicated that this is a site of an interaction between 
HB21 and a neighboring molecular cloud complex. In Figure [3] we present additional mosaicked 
ASCA GIS images which depict soft and hard emission (corresponding to the energy ranges of 
E<1 keV and E>1 keV, respectively) from this SNR. The GIS2 count rates for E<1 keV and E>1 
keV are 1.91(±0.08) x 10 -2 cts s _1 and 3.21(±0.15)xl0 -2 cts s _1 , respectively, while the respective 
count rates for E<2 keV and E>2 keV are 4.56(±0.13) x 1(T 2 cts s" 1 and 0.56(±0.08) x 10~ 2 cts 
s . Significantly more X-ray emission is detected from HB21 at E<2 keV X-ray energies than at 
E>2 keV energies, illustrating the soft spectral nature typical of SNRs. 

The GIS2/GIS3 spectra were extracted from elliptical regions approximately 23' in size care- 
fully selected to include most of the emission and avoid the edges of the field of view. Likewise, 
the SIS0/SIS1 spectra were extracted from square-shaped regions approximately 10 arcminutes on 
a side and again the edges of the fields of view were avoided. The spectral extraction regions 
are marked in Figure [0 We analyzed the extracted spectra using the software package XSP-EcQ 
Version 11.3.1 (|Arnaudlll996l ). For our spectral fitting we used two thermal models: the thermal 
model VAPEC, which describes an emission spectrum fr om a collisionally-ionized diffuse gas with 



variable elemental abundances (Smith et al. 



2000 



2001a 



b § and VNEI, which is a non-equilibrium 



collisional plasma model which assumes a constant temperature and single ionization parameter 
(jHamilton et al.l 119831 ; iLiedahl et al.1 Il995l ; iBorkowski et al.l l200ll ) . Photoelectric absorption along 
the line of sight was accounted for with the PHABS model; finally, we allowed the abundances of 
silicon and sulfur to vary during the fitting process (because lines associated with these particular 
elements are noticeable in the spectra) while leaving the abundances of the other elements frozen 
to solar values. 

In Table[2]we present results of our simultaneous fits to the GIS2/GIS3 and SIS0/SIS1 spectra 
for both the northwestern and southeastern regions of HB21. We have obtained statistically accept- 
able fits (with values of ~ 1.04-1.06) to the extracted spectra using a single thermal component 
(that is, either the VAPEC model or the VNEI model) for both the northwestern and southeastern 
regions. The column density and temperature are similar for both thermal models, namely Njj ~ 
2-3 xlO 21 cm" 2 and kT ~ 0.66-0.68 keV. In Figure 0J we present the extracted GIS2, GIS3, SIS0 
and SIS1 spectra for the northwest region of HB21: in each case, the fits obtained using ionization 
equilibrium (CIE) and nonequilibrium ionization models are comparable in quality. In addition, 
the abundances of silicon and sulfur in our spectral fits exceed solar abundances: these ASCA 
observations are the first to reveal enhanced a b unda nces of heavy elements in the X-ray spectra 
of HB21 (as also noticed by lLazendic &; Sland (|2006l )). Our results are consistent with previous 



4 http://heas arc. gsfc. nasa.gov. docs/xanadu/xspec/. 
5 Also see |http: / /hea-www.harvard.edu/ APEC 
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analyses of X-ray emission from HB21 ([Leahvlll987l ; lRholll995l ) where our present analysis of the 
broadband ASCA spectra has yielded similar temperatures to those derived from Einstein and 
ROSAT observations. However, with the data from the ASCA observations we may establish more 
stringent constraints on the ionization timescale and the abundances of sulfur and silicon. The ion- 
ization timescales derived with the VNEI model for the two regions are long (r = 5.9 (>3.2)xl0 11 



cm s and r 



4.1^;? xlO 11 cnT 



s for the northwestern and southeastern regions of HB21, re- 



spe ctively) and with i n the error bounds for this parameter, CIE is included (r > 10 cm s 



Smith fc Hughes! (|2010l )). For this reason and because both the VAPEC and the VNEI models 



see 

return equally acceptable fits, we argue that the X-ray emitting plasma located within the interior 
of HB 21 is close to ionization equilibrium. Similar to the results presented by lLazendic &: Slane 
(|2006l ). we also observed slightly enhanced silicon and sulfur abundances in HB21: based on our 
PHABSxVAPEC (PHABSxVNEI) fits, our measured abundances are Si=1.3±g:| (l-8±0.05) and 
S=2.4±1.0 (3.6ii;g) for the northwestern region and Si=1.4±0.3 (2.0±0.4) and S=1.7±g;| (3.0±1.4) 
for the southeastern region (see Table [2]). The slightly different abundances between the VAPEC 
and VNEI fits may be due to either differences in the predicted line strengths between the non- 
equilibrium and equilibrium conditions or different sets of atomic data incorporated in these models. 
The respective lower limits for the abundances of Si (S) for the northwestern region are 1.1 (1.3) 
and 1.4 (1.7), while the lower limits for the southeastern region are 1.1 (1.6) and 0.9 (1.6). The 
northwestern region shows modestly more consistent evidence of slightly enhanced silicon and sul- 
fur abundances than the southeastern region. This suggests a contribution to the observed X-ray 
spectra from ejecta material. 

In Figure [5] we present a plot of confidence contours for the silicon and sulfur abundances based 
on the PHABSxVAPEC fit to the spectrum of the northwestern region. No significant spectral 
differences are seen between the northwestern and southeas tern regions of this SNR. In their analysis 
of spectra extracted from the ASCA observations of HB21. lLazendic &: Sland (120061 ) found that fits 
obtained using the VNEI model to the spectra extracted for both regions gave better fits at a 
statistically-significant level (4cr) than fits obtained by a thermal plasma. In contrast, we find 
that fits obtai ned with VNEI and fits obtained with a standard thermal plasma are comparable in 
quality. Also, lLazendic Sz Sland (|2006l ) derived comparable (though slightly larger) values for Njj 
for both regions: those authors also found a similar trend where Njj is modestly elevated for the 
southeastern region compared to the north western region. We suspe ct that the minor differences 
between our results and those obtained by lLazendic fc Sland (J2006J) may be attributed to small 
differences in data re duction and spectral ana lysis techniques (such as background subtraction). 
We also point out that lLazendic &: Sland (|2006l ) commented that the derived column density values 
seemed to be rather high if HB21 is indeed only 0.8 kpc distant (which was their assumed distance 
to the SNR). They noted that an elevated column density (at least for the eastern portion of the 
SNR) is compa rable to the column den sity of the complex of molecular clouds seen toward HB21 
as described bv lTatematsu et al.l (|l990l ). We point out that if the larger distance to HB21 that we 
have assumed (1.7 kpc) is adopted, the measured values for N u seem more reasonable. 
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Lastly, to help determine more stringent constraints on the properties of the X-ray emitting 
plasma (specifically the abundances of silicon and sulfur), we simultaneously fit the GIS2 spectra ex- 
tracted for the northwestern and southeastern portions of the SNR with first the PHABSxVAPEC 
model and then with the PHABSxVNEI model. For both models we obtained statistically ac- 
ceptable fits (with x 2 f values of ~ 1.1) with values for kT and Njj that were consistent with 
those obtained for fits of the individual regions (that is, kT ~ 0.6 keV and Nh ~ 0.2 x 10 22 
cm -2 ). The abundances for both silicon and sulfur were indeed enhanced relative to solar: for the 
PHABSxVAPEC (PHABSxVNEI) model, the abundances were 1.7±0.4 (1.9±g;|) for silicon and 
(4.3^ 2 'g) for sulfur. The ionization timescale derived for the PHABSxVNEI model - r ~ 
4(>0.03) x 10 13 cm -3 s - is also consistent with ionization equilibrium. We summarize the results 
for these fits in Table [2j 



3.2. CTB 1 

A broadband (0.7-10.0 keV) exposure-corrected and mosaicked ASCA GIS image of CTB 1 is 
presented in both Figures [6] and [7] with radio contours from the CGPS overlaid. Both HB21 and 
CTB 1 show the typical center-filled X-ray morphology combined with a shell-like radio morphology 
that characterizes MM SNRs: in the case of CTB 1, however, as noted previously the X-ray emission 
is seen to extend through a gap along the northeastern portion of the radio shell. Figure [8] shows 
mosaicked images of the soft and hard X-ray emission (again corresponding to the energy ranges 
E<1 keV and E>1 keV, respectively); there is a noticeable difference in the spatial structure between 
the soft and the hard emission. In Figure [9] we present an optical Ha image of CTB 1 (courtesy of 
Robert Fesen) with the contours of the X-ray emission overlaid: the optical morphology of CTB 1 is 
quite simil ar to the radio mo rphology with the same shell-like structure and a prominent gap in the 



northeast (jFesen et al.lll997l ). Little Ha emission is seen in the interior of CTB 1, nor is any optical 
emission seen where the X-ray emission extends through the gap in the optical and radio shell in 
the northeast. The observed slight extension of X-ray emission in the west through the optical and 
radio shell is likely either residuals due to the broad ASCA point-spread-function or point sources 
rather than true emission from the SNR. The ASCA hard image reveals a point-like source seen in 
projection against the diffuse emission of CTB 1: it is located at RA (J2000.0) OO' 1 01 m 25. s 5, Dec 
(J2000.0) +62° 29' 40" and it lies very close to the eastern edge of the optical and radio shell. This 
source (denoted as 1WGA J0001. 4+6229) may be a neutron star possibly associated with CTB 1; 
we discuss it in detail in the next section. In contrast to HB21, a significant amount of emission 
from CTB 1 is seen at energies above 1 keV. 

Similar to our spectral analysis performed with HB21, we extracted GIS2 and GIS3 spectra 
from elliptical regions approximately 23' in size from both the southwestern and northeastern por- 
tions of the X-ray emitting plasma. When extracting GIS2 and GIS3 spectra from the northeastern 
region, we excluded a region 4 arcminutes in diameter centered on the position of the discrete X-ray 
source 1WGA J00001. 4+6229 to avoid spectral contamination by this source. We also extracted 
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SISO and SIS1 spectra from both the southwestern and northeastern portions of the SNR, again 
using square-shaped regions approximately 10 arcminutes on a side. Unfortunately, the signal-to- 
noise ratio of the extracted SISO and SIS1 spectra for the northeastern region was not sufficient for 
spectral analysis and we therefore omitted these spectra from our analysis. We attempted to fit the 
extracted ASCA spectra using the thermal models VAPEC and VNEI along with the model PHABS 
for the photoelectric absorption. In contrast to HB21, our fits with either VAPEC or VNEI did 
not produce a xt lower than 1.3 for either region and failed to account for the hard X-ray emission 
seen above ~3 keV from both regions. A two-component model with a soft thermal component 
with a temperature of 0.3-0.4 keV (either VAPEC or VNEI) and a second thermal component with 
a higher temperature or a power law component was required for an acceptable fit to the extracted 
spectra from both the southwest and northeast portions of the SNR (see Table U]). Thus, a hard 



excess is present in the ASCA spectra of CTB 1: lLazendic fe Sland (|2006l ) also identified a second 



thermal component of X-ray emission from this SNR based on analysis of extracted ASCA spec- 
tra but only for the southwestern region. This hard X-ray emission was not detected previously 
because the X-ray observatories employed in prior observations of CTB 1 (such as ROSAT) lacked 
the required sensitivity at high energies. 

We examine first the high-spatial resolution Chandra data before discussing the hard emission 
in more detail. In Figure [10] we present a three-color Chandra image of the interior X-ray emission 
surrounded by the well-defined optical and radio shell of CTB 1. To illustrate the spectral properties 
of this emission, in this Figure we have depicted soft (0.5-1.0 keV), medium (1.0-2.0 keV) and hard 
(2.0-8.0 keV) emission in red, green and blue, respectively. Many features with different spectral 
properties are visible: most importantly, features with primarily medium and hard spectra are 
clearly mixed together within the X-ray plasma. At the angular resolution of Chandra, the hard 
emission (such as the structure ~ 1 arcminute in size located at approximately RA (J2000.0) 23^ 
59 m 01.0 s , Dec (J2000.0) +62° 30' 57") is clearly diffuse and not point-like. 

To investigate variations in the spectral properties of the X-ray-emitting plasma of CTB 1 
as revealed by Chandra, we extracted spectra from three different Chandra chips. These three 
extraction regions are as follows: the first region (which we refer to as the "diffuse" region) is on 
the ACIS-I2 chip and covers most of the area of this chip. A second region (which we refer to as the 
"soft" region) is located on the ACIS-I3 chip and covers most of the area of that chip as well. Lastly 
we extracted spectra from a third region (which we will call the "hard" region) which corresponds 
to a region of hard emission mentioned above and is located on the ACIS-I1 chip: the positions 
of all three regions are indicated in Figure [TUl A region of excess medium energy seen toward the 
middle of the field in Figure [TUl largely fell into gaps between the Chandra detector chips: for this 
reason, a detailed analysis of its X-ray spectrum could not be conducted. 

We present the extracted spectra of all three regions in Figure [TTJ spectral variations are 
present across the X-ray emitting plasma of CTB 1. For example, the spectrum from the "diffuse" 
region shows prominent O Ly a (0.65 keV), Ne IX (0.9 keV), and Mg XIII (1.35 keV) lines, 
together with Fe L-shell line emission: in this spectrum the O and Ne lines are stronger than 
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the Fe L-shell line emission. In contrast, the lines are hardly noticeable in the "hard" region 
spectrum. We derived an acceptable fit to the spectrum of the diffuse region using a thermal model 
(VAPEC) with a temperature kT = 0.28±0.03 keV and N H = 0.64±0.08x 10 22 cm" 2 ; this fit reveals 
enhanced oxygen and neon abundances and a low iron abundance (see Table [3]). The column density 
-/Vh-=0.64± .08 x 10 22 cm -2 is consistent with the optical extinction of E(B-V) =0.7-1 derived by 



Fesen et al.1 (j 19971 ). which is equivalent to iY#~0. 5-0.7 x 10 cm . In Figures fl~2a and [T2b. we 



present confidence contour plots for the abundances of oxygen compared to iron and neon compared 
to iron (respectively) based on the fit derived from the PHABSx VAPEC model to the spectrum 
of the "diffuse" region. These figures show that the oxygen and neon abundances are both above 
solar while the iron a bundance is below solar: such relative abundances are a typical characteristic 
of oxygen-rich SNRs (jWooslev fe Weaver! 1 19951 ). 



Our detection and analysis presented here of probable oxygen-rich ejecta in CTB 1 is the 
first detailed stud y presented of oxygen-ri ch ejecta associated with an MM SNR. Based on X-ray 



spectral analysis, lLazendic fe Sland (120061 ) suggested that another MM SNR - HB 3 - may also 



feature oxygen-rich ejecta, though those authors could not determine if the plasma associated with 
that SNR has significantly enhanced abundances of oxygen, neon and magnesium or marginally 
enhanced abundances of magnesium and underabundant iron. The Chandra spectra are not well 
fit by a single temperature thermal model (see the poor match to the Mg XII Lya line as seen in 
Figure [3]), so it is not possible to have oxygen, neon and magnesium all residing in the same constant 
temperature plasma in ionization equilibrium. Magnesium appears to be more ionized than neon. 
This may be caused by higher temperatures in magnesium-rich ejecta than in the oxygen- and 
neon-rich ejecta. It is also possible that the high temperature plasma could be ove rionized; an 
overionized plasma has been reported in the MM SNR IC 443 (IKawasaki et al.ll2002l ). When we 
fit the spectrum of the "soft" region using a VAPEC model combined with the PHABS model for 
the interstellar photoelectric absorption, the iron abundance is higher and the neon abundance is 
lower when compared with the "diffuse" region spectra (see Table [3j), suggesting spatial variations 
in the chemical composition within the X-ray emitting gas of CTB 1. 

Lastly, we describe the spectral analysis of the third ("hard") region which features a harder 
spectrum when compared to the other two spectra (see Figure fTTj) . Such localized hard emission 
has been reported in several oth er MM SNRs, either as n on-thermal diffuse knots seen in XMM- 
Newt on observation s of IC 443 (IBocchino fe Bykovi 120031 ) or as a pulsar wind nebula (also in IC 



443 - lOlbert et al. I (|200ll )). Because of the lack of lines in the spectrum, we first fit the spectra 



with a power law, which resulted in an acceptable fit, but with an unrealistically high photon index 
of r~6.4l|j (see Table©. In comparison, a VAPEC thermal model yielded a comparable-quality 
fit with a temperature of kT = 0.66^q' 2 J keV: for this fit, we froze the abundance of oxygen to 
1.7 (to be consistent with the fits derived to the spectra of the "diffuse" and "soft" regions) or 1 
while fitting for the abundances of neon and iron. We also fit the spectrum of this region with a 
combination of an APEC thermal component and a non-thermal power-law component; to reduce 
the number of free parameters we fixed the temperature of the thermal component to kT = 0.28 
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keV, equal to the temperature derived from fits to the other two regions. The results are given 
in Table [3j unfortunately, we do not have enough counts in the spectrum of the "hard" region 
to distinguish between non-thermal and thermal interpretations of the spectrum of this source. 
However, the temperature derived from the thermal fit (kT = 0.66 keV) is significantly higher 
than the temperature of the "soft" and "diffuse" regions. Also, thermal models systematically 
underpredict Chandra spectra at high energies in all the regions: we interpret this result as evidence 
that the hard excess observed toward CTB 1 has a diffuse origin. 

In Figure [T3l we present a plot of confidence contours for N h arid kT that correspond to 
the fit to the spectrum of the "hard" region as fit with the PHABSxVAPEC model with fixed 
solar abundances. The spectrum of this region demonstrates a bimodality with temperatures of 
kT ~ 0.28 and kT ~ 0.66 keV, suggesting the presence of an additional component in addition 
to the thermal component associated with the "diffuse" region. We will discuss the nature of this 
hard component further below and in Section [7J We also note diff erences between our resu lts from 



analyzing Chandra ACIS-I spectra and the results presented by lLazendic fe Sland (|2006l ): those 
authors jointly fit four spectra taken from each ACIS-I chip (they did not attempt any spectral 
analysis of small regions like the "hard" region") and presented an acceptable fit obtained using 
two VNEI components: one with a temperature kT = 0.201q'q^ keV and solar abundances and 
the other with a temperature A:T=0.86to'og and an elevated magnesium abundance (Mg=3.li ; 4 ). 
Those authors did comment on the presence of lines associated with oxygen as well as the neon and 
iron line blend in the extracted ACIS spectra but they did not present an analysis of the abundance 
of those elements. 

After establishing the spectral properties of the soft X-ray emitting plasma from Chandra 
spectra, we fit the ASCA GIS2/3 and SIS0/1 spectra of the southwestern portion of CTB 1. 
We froze the oxygen abundance to 1.7 (or 1) for all spectral fits while allowing the neon and 
iron abundances to vary. Because models with a single thermal component were not sufficient to 
describe the X-ray spectra, we used a combination of the thermal models VAPEC and VNEI along 
with a power law model to jointly fit the ASCA spectra: the results of these spectral analyses are 
given in Table [H From these fits, we estimate a column density Nh ~ 0.5-0.6 x 10 22 cm" 2 and a 
temperature kT so { t ~ 0.2-0.3 keV for the soft emission. The soft thermal component has most likely 
attained thermal equilibrium because fits to this soft emission using the VNEI model resulted in 
a long ionization timescale (r ~ 1 (>0.2) x 10 11 cm -3 s). The hydrogen column density Nh and 

"anent kT so f t are in agreement with the analysis of R OS A T spectra 



(Rho 


1995; 


Craie et al. 


1997) 



statistically-acceptable fits to the ASCA spectra of the southwestern region (with values for the xt 
~1.1): the addition of a power law with a photon index T ~ 2-3 or a second thermal component 
with a temperature kT ~ 3 keV yields fits with a comparable quality. Unfortunately due to a low 
number of counts in the spectra at the higher X-ray energies, we cannot distinguish between different 
models for the high-energy emission. In Figure [14] we present the GIS and SIS X-ray spectra of 
the southwest region, fit with the combination of a thermal (VAPEC) and non-thermal (power 



-14- 



law) model. Our results differ from lLazendic fe Sland (|2006l ) who fit the extracted ASCA spectra 



from this region with two thermal components in CIE: the softer temperature component featured 
a temperature kT = 0.19^o'q3 keV and a magnesium abundance fixed at solar while the harder 
temperature component featured a temperature fcT=0.82^QQg keV and an elevated magnesium 
abundance (Mg=2.7±£§). 

Finally, we examine the ASCA spectra of the northeast region of CTB 1 which corresponds to 
the known "break-out" site seen in optical and radio images of this SNR.. Like the southwest region, 
a second component is needed (in addition to a soft thermal component) to derive a statistically 
acceptable fit. A fit to these spectra using a power-law component for the hard emission is presented 
in Table 2] and the GIS spectra are presented in Figure [15j For the thermal component we have 
first assumed the abundances of oxygen, neon and iron to be 1.7, 1.6 and 0.4, respectively, equal 
to the abundances in the "diffuse" region. Secondly, we assumed solar abundances for the hard 
component. Although the GIS spectra of the northeastern region feature a stronger Fe L-shell line 
complex when compared with the spectra of the southwestern region, our spectral fits could not 
confirm non-solar abundances and the abundances are consistent with solar values. 

The photon index derived from fits to the northeast is flatter than the photon index derived 
from fits for the southwest region: T = 1.4 compared to T = 2 — 3, respectively. However, the 
photon statistics at higher energies is poor, making it difficult to determine the true nature of this 
hard emission. Although the northeast region does correspond to the prominent breakout feature, 
we do not find evidence for any major differences in the X-ray properties between the northeast and 
southwestern regions of CTB 1. The fact that no significant variations are seen on large scales in 
the X-ray properties of CTB 1, coupled with the variations seen on small scales as revealed by the 
Chandra observation, indicate that the X-ray emission from this SNR is complex. A parallel may 
be drawn with the results from Chandra observations of 3C 391, another MM SNR: in the case of 



that source, local spectral differences appeared to be stronger than global ones (IChen et al.ll2004f ). 



Here again we find that our results differ from those presented by lLazendic fe Sland (|2006l ) : for 
this region, the authors derived an a 
keV) in CIE with solar abundances. 



this region, the authors derived an acceptable fit using a single thermal component (kT = 0A8^q'q® 



4. 1WGA J0001.4+6229 An X-ray Pulsar Associated with CTB 1? 



The ASCA hard energy image (E > 1 keV) revealed a hard source in the northeastern region 
of CTB 1 which is located just inside the eastern shell of the SNR: the position of this source is 
RA (J2000.0) 00^ 01 m 25. s 5, Dec (J2000.0) +62° 29' 40" with a positional uncertainty of 13". This 
i s the discrete X-ray sou r ce 1 WGA J0001. 4+ 6229 in the Catalog of ROSAT PSPC WGA Sources 
(|white et al.1 Il994l . Il997l : kngelini et all boocl fl This X-ray source may possibly be a neutron star 



6 Also see 



http: / /wgacat. gsfc.nasa.gov 
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associated with CTB 1. We therefore conducted spectral and timing analysis of the X-ray emission 
from this source using the ASCA GIS2 and GIS3 datasets. 

The procedure for extracting GIS2 and GIS3 spectra of 1WGA J0001. 4+6229 and performing a 
spectral analysis was the same as for extracting GIS2 and GIS3 spectra of the diffuse emission from 
HB21 and CTB 1: a circular region four arcminutes in diameter centered on the source position 
was used to extract spectra. The total number of counts and the corresponding count rate (over 
the energy range from 0.6 keV to 10 keV) for our GIS2 and GIS3 observations were 115 and 103 
counts, and 5.71+0.53 x 10~ 3 and 5.11+0.51 xlO -3 counts per second, respectively. We derived a 
statistically-acceptable joint fit to the spectra using a simple power law model combined with the 
same PHABS model mentioned previously for the photoelectric absorption along the line of sight. 
The parameters of this fit were a column density of Nh=0.3 (<0.65) x 10 22 cm -2 and a photon 
index r=2.2l ( j ) '2 : m Figure [TBI we present the extracted GIS2/GIS3 spectra together with the best- 
fit model, and a confidence contour plot for Njj and T. The derived photon index is typical for 
rotation-powered pulsars. The column density is consistent with the range of column densities 
derived in our fits (see Tables [3] and H]) to the CTB 1 spectra, hinting at a possible association. 
If we fit the extracted spectra for this source with a blackbody model, a temperature of kT~l.l 
keV is derived (although this fit with x 2 =0.92 for 48 degrees of freedom (DOF) is inferior to the 
power-law fit with x^=0.79). Because of the low derived column density, our estimated absorbed 
and unabsorbed fluxes for this source are virtually identical: for the GIS2 (GIS3) spectrum, the 
flux is 5.4xl0 -13 (6.4xl0 -13 ) ergs cm -2 sec -1 ; at the assumed distance to CTB 1, these fluxes 
correspond to luminosities of 6.2xl0 32 (7.4xl0 32 ) ergs s _1 , respectively. We also performed a 
timing analysis using the GIS2/GIS3 datasets to search for pulsed X-ray emission from this source 
and detected a period of 47.6154 milliseconds using the Rayleigh test (a maximum signal of Z 2 = 
31.4), but the detection is not statistically significant. 

We further searched for pulsations from 1WGA J0001. 4+6229 with the 100- meter Green Bank 
Telescope (GBT) of the National Radio Astronomy Observatory 

(NRAO0) 

on 2004 December 7. 

The target position was observed for 7.6 hours at a center frequency of 825 M Hz. The frontend was 
the GBT Prime Focus 1 receiver to feed the Pulsar Spigot ( Kaplan et alJ[2005l ) and Berkeley-Caltech 



Pulsar Machine (BCPM) backends. The receiver provided 50 MHz of bandwidth in two orthogonal 
polarizations that were summed and synthesized into 1024 frequency channels every 81.92 (is in the 
Spigot and 96 channels of 0.25 MHz width every 144 //s in the BCPM. The interstellar dispersion 
toward CTB1 is unkno wn, but we can estimate it with the latest model of Galactic electron density 



( Cordes fc Lazid 120021 ) . which predicts a dispersion measure (DM) of 105 pc cm -3 for a distance 
of 3.1 kpc or DM = 33 pc cm -3 at a distance of 1.6 kpc (corresponding to the two distances to 
CTB 1 that have been published in the literature). We therefore take a conservative upper limit of 
DM = 1000 to CTB 1 (given the narrow channels and relative long pulsation period the search is 



7 The National Radio Astronomy Observatory is a facility of the National Science Foundation, operated under 
cooperative agreement by Associated Universities, Inc. 
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not highly DM dependent). The data set was dedispersed with DMs from to 1000 and searched 
for periodicities using standard folding and fast Fourier- Transform (FFT)-based techniques. Based 
on these analyses, we find no significant evidence for pulsations with any period from 1WGA 
J0001.4+6229. 

Assuming that 1WGA J0001. 4+6229 is in fact a neutron star associated with CTB 1, a trans- 
verse velocity can be estimated. The angular displacement of 1WGA J0001. 4+6229 from the center 
of CTB 1 is 14' while the radius of the SNR itself is 17'. Therefore, we calculate a transverse velocity 
v = 850 d .s.i t] \ km s" 1 , wh ere d3.i is the distance to CTB 1 in units of 3.1 kpc and ti.6 is the age of 



the SNR ( Fesen et al,lll997l ) in units of 1.6 x 10 4 yr. This estimated transverse velocity is high but 
this may be an overestimate because of the considerable uncertainties associated with estimates of 
the distance and age of CTB 1: if we assume a distance to the SNR of 1.6 kpc, the tranverse velocity 
is only 420 km s _1 . In particular, the published age estimates are based on simple one-dimensional 
SNR models: the obvious breakout morphology of CTB 1 clearly indicates that such models are 
not applicable in this case. For comparison, the transverse velocities for ne utron stars locate d off 
center in their associated SNRs are 375 km s _1 i n the case of the SN R W44 ( Frail et al.lll996l ) and 



250+50 km s" 1 in the case of the SNR IC 443 (IQlbert et al. 



2001 



It is plausible that 1WGA 
J0001. 4+6229 is associated with CTB 1 but deeper high-spatial resolution X-ray observations are 
needed to examine its spectrum in more detail and search for possible pulsations. 



5. Plasma Conditions in HB21 



We first estimate the density and mass of the X-ray emitting plasma associated with HB21 
based on the emission measures derived from our spectral fitting. Our GIS spectral extraction 
regions extend over approximately 11'5 x 11^5 or 5.7 pc x 5.7 pc (2.7 pc x 2.7 pc) at the assumed 
1.7 (0.8) kpc distance to HB21. Assuming a cylindrical geometry with the long axis equal to the 
observed extent of the X-ray plasma (35'8, corresponding to 17.7 (8.3) pc), the volume of each 
region is approximately 5.4xl0 58 (5.6 xlO 57 ) cm 3 . From the mean values of our derived emission 
measures (which are approximately the same for all regions and all models) , we calculate an electron 
density n e ~ 0.06 (0.08) cm -3 (where we have assumed n e ~ 1.2ra#) and a volume filling fraction 
of unity based on the smooth appearance and isothermal nature of the X-ray emitting gas. Based 
on this value, we estimate the total mass of the X-ray emitting plasma within the field of view 
of the ASCA observations to be only ~2.6(1.5)Mq. When we account for the incomplete spatial 
coverage of HB21 by the ASCA observations, the total X-ray mass could be higher by a factor of 9, 
which amounts to a total of 23.4 (14) M@. The corresponding Si and S masses are 1.7 x 1O" 3 M 
and 1.3 x 10~ 3 Mq, respectively. The presence of a bright radio shell without associated X-ray 
emission combined with t he detections of an expanding HI shell and in frared emission from shock- 



cloud interaction regions (|Koo h Heiles 



1991 



Shinn et al.l |2009i . l2010l . and the references therein) 



imply that HB21 is in a radiative cooling stage. The SN R age inferred from the presence of an 



expanding 



J exp 



124 km s" 1 ) HI shell is t d = 4.5xl0 4 yr (jKoo fc Heiles 



1991 



We can also infer 
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the pressure within HB21 from properties of the X-ray emitting gas. The total number of particles 
is n tota [ = n e +n#+n# e ~2n e for a plasma with cosmic abundances: from our estimated values for 
electron density and temperature, the corresponding pressure is P/k = 2n e T = 0.9 (1.2) xlO 6 K 
cm -3 , which as about two orders of magnitude higher than the typical ISM pressure. The physical 
properties of HB21 are summarized in Table 

The X-ray properties of HB21 are similar to those measured for many other MM SNRs. First, 
the presence of an isothermal plas ma with a tempera t ure of kT ~ 0.2-0.7 keV is consistent with 



other MM SNRs s uch as 3 C391 dRho fc Petrel Il996l : Ichen et alJboO^h . W44 faho et al 



20041 ). and Kes 79 (|Sun et al 



Shelton et al.ll2Q04h 3C400.2 (Yoshita et all200lh . W51C (jKoo et alJ2002h . W63 (jMavromatakis et al 



1994 



2004! ) . This result supports the interpretation that these SNRs are 



evolved and in the radiative phase as suggested by the presence of infrared, optical and HI shells 
for many of these SNRs. There are no temperature variations and no pronounced enhancements 
of chemical abundances in HB21 either, just as in many other MM SNRs like 3C 391. These 
properties of HB21 exemplify the typical X-ray properties characteristic of MM SNRs as defined by 
Rho fc Petrel ( 19981 ). At a sufficiently old age (~ 10 6 yr) age, a S NR should exhibit a centrally-filled 
X-ray morphology and eventually merge with the hot ISM gas (jCui &; Cox! 1 19921 ): however, MM 
SNRs attain this state at a much earlier age (~ 10 4 yr). When a distance of 1.7 kpc to HB21 
is assumed, the calculated X-ray emitting mass of this SNR is comparable to those of other MM 
SNRs. For standard radiative SNR models, we exp ect ~ 100 Mp of X-ray emitting gas at the 



HB21 age of 4 x 10 4 yr ba s ed on equations given by ICui Coxl ( 19921 ) and models presented by 
Hellsten fc Sommer-Larsenl (119951 ). Even more X-ray emitting gas is expected in conduction models 



Cox et al.1 (|l99flh . This discrepancy betwee n the observed and predicted X-ray emitting mass is 



of 

also present in W28 (IRho fc Borkowskil 120021 ) . but unlike in W28 (where large temperature gradi 
ents have been detected), the presence of an is othermal plasm a at the center of HB21 su ggests that 



the ele ctron thermal conduction is important (jChevalier 



1999). The conduction model of lCox et al 



(| 19991 ) overpredicts the mass of X-ray emitting gas: however, this model assumes a uniform ambient 



ISM while HB21 i s kno wn to be interacting with with clumpy molecular clouds (jKoo et al.l 12001 



Shinn et al.ll2009i . l2010l ). More elaborate X-ray emission models of SNRs in molecular clouds are 
needed to account for the observed X-ray properties of HB21 and similar MM SNRs. 



6. Supernova Ejecta in CTB 1 



The enhanced abundances of oxygen and neon and low iron abundances in the "diffuse" region 
(see Table [3]) indicate that CTB 1 is likely an oxygen-rich SNR. This SNR was likely produced by 
a core-collapse SN explosion, because such exp losions produce O- and Ne-rich, and Fe-poor ejecta 
( Nomoto et al.lll997l : IWooslev fe Weaver! Il995l ). This finding is consist ent with the presence of a 



massive star forming environment near CTB 1 (jLandecker et al.lll982l ): in addition, the scenario 
for the creation of this SNR by a core-collapse SN explosion would be further supported if the 
discrete X-ray source discussed earlier - 1WGA J0001. 4+6229 - is shown to be a neutron star as- 
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sociated with this SNR. Examples of oxygen-rich SNRs include young sources like Cas A, N132D, 
and E0102.2 -72.3; recently two older SNRs located in the Small Ma gellanic Cloud (S MC) - SNR 
B0049-73.6 (|Hendrick. Reynolds, k Borkowskilhood ) and B0103-72.6 (|Park et allbood ) - were also 
classified as oxygen-rich SNRs. Both of these SMC SNRs show the ejecta material in their interiors 
surrounded by shells of swept-up ambient material at relatively low X-ray emitting temperatures. 
An X-ray emitting shell might be present in CTB 1, but its detection may be prevented by sub- 
stantial interstellar absorption in this direction combined with an expected low temperature of the 
shocked ambient gas. 

We estimated the X-ray mass and density of CTB 1 from our fits to the ASCA spectra, 
assuming metal abundances derived from the Chandra spectra of the "diffuse" and "soft" regions. 
The spectral fits to the GIS spectra in the ll'.h x llf5 region imply an electron density of 0.16/ so y f 
cm -3 for the soft thermal component (where f so ft is the volume filling factor for this component). 



-1/2 



The corresponding hydrogen density nn is equal to n e /1.2 = 0.13/ so j t 
we estimate the total mass of the X-ray emitting plasma to be ~ 40J so j t 



cm 3 ; based on this value 



1/2 

40/ 'r,M Q . From our derived 



abundances of oxygen, neon and iron based on Chandra data, we estimate the oxygen, neon and 



iron masses to be 0.66/^^, 0.11/^^, and 0.03/^^ M , respectively. The ratio of [O/Fe] is 4.3+2 1 
and [Ne/Fe] is 4.0l|;° for CTB 1. The expected ratio of [O/Fe] is 0.75 for a Type la explosion and 
greater than 4 for a core-collapse explosion. These abundances imply that CTB 1 is a remnant of a 
core-collapse e xplosion and are consistent with the prediction s for a stellar progenitor with a mass 
of 13 - 15 M Q (jWooslev &: Weaver! Il995r , iNomoto et al.lll997l ). but higher mass stellar progenitors 
are not excluded. 

Finally, we estimate the pressures of the soft and hard components of the X-ray emitting gas 
using the parameters of the PHABSx(VAPEC+VAPEC) model: for the soft component we calcu- 
late a corresponding pressure P/k = 1.1 xlO 6 f S a(t K cm_3 - F° r the hard component (assuming 
a thermal origin), we first need to calculate the corresponding electron density which we can de- 
termine from the electron density of the soft component and the emission measures (EMs) of the 
soft and hard components (i.e., n e (hard) = ra e (soft) [EM/j^d/EM^o^] 1 / 2 ). From this relation, we 
obtain n e (hard) = 0.029 //J/J cm_3 (here fhard is the volume filling factor for the hard component) 
and therefore a corresponding pressure P/k = 2.0xl0 6 fullj K cm -3 . This result implies a factor 
of four larger filling factor for the hotter gas than the cooler gas if these two components are in 
pressure equilibrium: a higher filling factor for the hot gas is typical. Assuming f so ft + fhard = 1 5 
pressure within CTB 1 is 2.8 x 10 6 K cm~ 3 . We summarize these inferred physical properties for 
CTB 1 in Table [5j 

CTB 1 therefore belongs to a growing number of known evolved SNRs which feature an 
enhanced metal abundance in their interiors. An exam ple of another MM SNR which features 
such enhanced abundan ces is W44 (jShelton et al.1 [2004]) : other similar sources are identified by 
Lazendic fc Sland (|2006l ) (including HB21, which was analyzed both in their stud y and in the 



,1/2 



1/2 



j+10.2 



study presented here.) In addition, two other Galactic SNRs - the Cygnus Loop (jMivata et al 
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199a ) a nd G347.7+0.2 (lLazendic et al.l 1200a ) - feature enhanced abundances of metals as well. 
W49B rtHwang et al.l 120001 ) sh ows highly enhance d abundances but its age is estimated to be 2000 
years ([Hwang et al.ll2000l ) and lRho &; Petrel (|1998l ) describe the source as an atypical MM SNR. Be- 
cause MM SNRs lik e CTB 1 are comm only believed t o be evolved sources - age estimates of CTB 1 
range from 9000 yr ( Craig et al.lll997) to 4.4xl0 4 vr (IKoo &: Heileslll99ll ) - their X-ray spectra are 
dominated by swept-up material (|Rho &: Petrel Il998l ). Therefore, the detection of X-ray-emitting 
material associated with these sources with enhanced metal abundances is unexpected. The de- 
tection of O-rich ejecta associated with CTB 1 is particularly noteworthy: CTB 1 may belong to 
a previously unrecognized class of MM SNRs whose X-ray emission is dominated by O-rich ejecta 
located within their interiors. As noted previousl y, another possible mem ber of this particular 
class of MM SNRs with O-rich ejecta may be HB3 (jLazendic fc Sland 12000 ) . We note that O-rich 
ejecta has been previously detected in the evolved (1.4 xlO 4 yr old) SMC SNR B0049-73.6 by 
Hendrick. Reynolds. fe Borkowskil (|2005l ) and it is likely that centrally-located ejecta will be found 
in a number of relatively old Galactic SNRs. 

The X-ray emitting pla sma associated with C TB 1 clearly ex tends through the gap in the 
crescent-shaped radio shell: lHailev &: Craig (|1994l ) and iRhol (|1995l ) first noticed this remarkable 
extension of X-ray emission based on ROSAT PSPC observatio ns. Two scenarios have been pro- 
posed to explain the morphology of the observed X-ray emission: ICraig et al.l (j 19971 ) has suggested 
that the ambient ISM toward the northeastern portion of the SNR was cleared by a supernova event 
which took place prior to the birth of CTB 1, and thus a breakout occurred as the SNR expanded 
into this reg ion of a dramatically lowe r density. A competing theory for the morphology has been 
proposed by lYar-Uyaniker et al.l ([2004]), who suggested that the X-ray emission from CTB 1 lies in 
the interior of a bubble seen in the 21 cm H line, presumably blown by winds of the CTB 1 stellar 
progenitor. Our X-ray images show that the diffuse X-ray emission in the northeast clearly extends 
through the relatively narrow break in the optical and radio shell, with the breakout directed into 
the interior of the bubble seen in the neutral hydrogen line. Such a morphology favors a scenario 
where a supernova explosion occurred within the HI shell, followed by subsequent breakout into 
the bubble and not an explosion within the bubble itself. Our ASCA spectra show little difference 
between the X-ray properties of the southwest and northeast regions: hints of variations in temper- 
ature and abundances exist but better X-ray data are needed to determine whether they are real 
and not just statistical fluctuations. 



7. The Nature of the Hard X-ray Emission from CTB 1 

The ASCA spectra of the southwestern portion of CTB 1 revealed the presence of a hard 
component in addition to the soft (kT ~ 0.28 keV) component. This hard component may be 
modeled as a second thermal component with a temperature of kT ~ 3 keV or as a power-law 
continuum with a photon index T ~ 2-3. Hard X-ray emission was also detected by ASCA in the 
northeast region of CTB 1: a power-law component with a photon index r ~ 1.4 (a somewhat 
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lower value compared to the southwest region) combined with a soft thermal component (again 
with a temperature kT ~ 0.28 keV) yields a statisically acceptable fit. The Chandra observation 
of the southwest region of CTB 1 revealed regions of harder emission patches on the scale of an 
arcminute in size: the spectrum of one of these regions can be modeled by either a single thermal 
component with an elevated temperature (kT ~ 0.66 keV) or as the combination of a soft thermal 
component (kT ~ 0.28 keV) and a power law component with a photon index T ~ 2.0 (see Tables 
[3] and [4]). The "hard" region observed by Chandra is only a few arcminutes from the center of CTB 
1 and well inside the radio-emitting shell of the SNR: it is diffuse in nature although the number of 
counts detected from the source is limited. It is possible that the hard X-ray emission detected by 
ASCA from CTB 1 may be composed of localized hard regions such as this one: unfortunately we 
do not have enough counts in this "hard" region to distinguish between thermal and non-thermal 
origins. We note that two other MM SNRs, W28 and IC 443, contain high-temperature thermal 



plasmas in their interiors (jRho &; Borkowski 



2002 



Kawasaki et al.ll2002l ). 



Several possible explanations may be considered for the origin of hard X-ray emission from 
MM SNRs: first, the hard emission may be caused by temperature variations within the SNR. 
In the case of CTB 1, this scenario is supported by a good fit to the spectrum of the "hard" re- 
gion with a thermal component with a much higher than average temperature. Supernova ejecta 
may be inhomogeneous, in which case a multi-temperature plasma with spatially-varying abun- 
dances is expected. Alternatively, the "hard" regions may be caused by localized non thermal emis- 
sion: such emission ha s already been detected in IC 443 (jBocchino fc Bvkovl |2003| ) and 7 Cygni 
(jUchiyama et al.ll2002i ). Additional observations are needed to understand the true nature of the 
hard X-ray emission from CTB 1. 



8. Summary 

1. We presented ASCA observations of the MM SNR HB21. Our ASCA images of this SNR are 
similar to R OS A T images and reveal a diffuse centrally filled X-ray emission located within a radio 
shell. Prom X-ray spectra, we measure a column density toward this source and a temperature for 
the X-ray emitting plasma of N h ~ 0.3 xlO 22 cm -2 and kT ~ 0.7 keV, respectively: no significant 
spatial differences in temperature are found. Silicon and sulfur abundances are slightly enhanced 
relative to solar, particularly for the northwestern region, and no hard component to the X-ray 
emission was detected. The properties of HB21 are similar to those seen in several other MM 
SNRs, such as the presence of isothermal plasma. This result supports the interpretation that MM 
SNRs are evolved sources currently in the radiative phase of evolutio n: the X-ray proper ties of 



HB21 exemplify the primary characteristics of MM SNRs as defined by lRho Sz Petrel (| 19981 ). 



2. We presented ASCA and Chandra observations of the MM SNR CTB 1. ASCA observations 
reveal center-filled X-ray emission located within the radio shell: the X-ray emission extends outside 
the circular shell through the breakout gap in the northeast. While the global X-ray and radio 
morphology is similar to HB21, the X-ray spectra of CTB 1 and HB21 are very different. The 
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X-ray spectrum of CTB 1 shows several prominent lines such as O Lya (0.65 keV) and Ne IX (0.9 
keV). We find that CTB 1 is likely an oxygen-rich SNR with enhanced abundances of oxygen and 
neon: this is surprising for an evolved SNR such as CTB 1. The derived abundances are consistent 
with an explosion of a stellar progenitor with a mass of 13 - 15 M Q and possibly even higher. 

3. The ASCA spectra of the southwest region of CTB 1 cannot be fit with a single thermal 
component and instead require the presence of an additional component to account for an excess 
emission seen at higher energies. Based on ASCA and Chandra spectra of CTB 1, we derive a 
column density N h ~ 0.6 x 10 22 cm -2 and the soft component temperature feT so /t~ 0.28 keV; the 
hard emission may be modeled either by a thermal component with a temperature kT^ard ~ 3 keV 
or by a power law component with a photon index of T ~ 2-3. Likewise, the ASCA spectra of the 
northeast region of CTB 1 also show an excess at higher energies: these spectra are fit best by a 
power law with a photon index T = 1.4 plus the soft thermal component. The Chandra observation 
of the southwestern region reveals localized regions of hard emission: one such region is ~ 1' in size. 
The X-ray spectrum of this region may be fit with either a higher temperature thermal component 
(kT = 0.66 keV) or with the combination of a softer thermal component (kT = 0.28 keV) and a 
power law component (F ~ 2.0). Because of the poor photon statistics, its true nature is unclear. 
Possible scenarios for its origin include temperature variations within the X-ray emitting plasma 
of CTB 1, including the ejecta, or localized non-thermal X-ray emission. 

4. The ASCA hard (E > 1 keV) image of CTB 1 reveals a point-like source seen in projection 
against the diffuse emission of CTB 1. This source - denoted as 1WGA J0001. 4+6229 and located 
at RA (J2000.0) 00 h 01 m 25. s 5, Dec (J2000.0) +62° 29' 40" - may be a neutron star associated 
with CTB 1. The GIS2/GIS3 spectra of this source are well-fit by a power-law continuum with a 
photon index r=2.2^ < j > 2 (typical for rotation-powered pulsars) and the measured column density is 
comparable to the column density measured for CTB 1. There is marginal evidence for pulsations 
in X-ray data at 47.6 msec, but no pulsations have been detected at radio wavelengths. 
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Table 1. Summary of ASCA GIS and SIS Observations of HB21 and CTB 1 
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Note. — The units of Right Ascension are hours, minutes and seconds and the units of Declination are degrees, arcminutes and arcseconds. Count rates are for the energy range 0.6-10.0 
keV. 
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Table 2. Summary of Fits to GIS and SIS Spectra of HB21 a 



Parameter 


GIS2/3+SIS0/1 
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GIS2 


GIS2 


Region 


Northwest 


Northwest 


Southeast 


Southeast 


Northwest 
and Southeast 


Northwest 
and Southeast 


Model b 


PHABSx 


PHABSx 


PHABSx 


PHABSx 


PHABSx 


PHABSx 




VAPEU 


VNEI 


VAPEC 


VNEI 
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a All quoted error bounds are 90% confidence intervals. 

b PHABS is a photoelectric absorption model, VAPEC is a thermal plasma model in ionization equilibrium, and VNEI is a noncquilibrium ionization thermal 
model (see Section |3. II for references for these models). 

c Emission measure, defined here as Jn e ri[f dV: here, d is the distance to HB21 (in cm) and n e and riu are the electron and H densities (in cm -3 ). 

d For the energy range 0.6-10.0 keV. The luminosity estimates are for an assumed distance of 1.7 kpc. 



Table 3. Summary of Fits to Chandra Spectra of CTB 1 
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Hard Region 


Model 


PHABSx 


PHABSx 


PHABSx 


PHABSx 


PHABSx 




VAPEC 


VAPEC 


VAPEC 


POWER LAW 


(APEC+POWER LAW) 


Xl (x 2 /DOF) 

N H (x 10 22 cm- 2 ) 

fcTi (keV) 

O 

Ne 

Fe 

EMi b /(47rd 2 /l(r 14 ) (cm- 5 ) 

r 

Normalization 


0.99 (169.72/172) 
0.64±0.08 
0.28±0.03 

1 -'-0.6 
J - D -0.B 

0.4±0.2 
4.33X10- 3 


0.94 (122.77/131) 
56+ 011 

n no + 0.06 
°' 28 -0.04 

1 s+ 2 - 2 

L - L -0A 
n 7 +0.6 

2.23x10-3 


0.41 d [0.41] e 
0.18«0.75) [0.36lE!'39] 
0-66+ 41 [0.63+ 53 ] 
1.7 [=1] 

i-i±S:2 m 

0.5±2;° [=1] 
2.14X10- 4 


0.50 d 
0.76(1°^) 

2.08X10" 4 


0.40 d 
0.47(-) 
0.28 (frozen) 
1 (frozen) 
1 (frozen) 
1 (frozen) 
2.41X10- 4 

1.9(-) 
1.27X10" 5 


Counts 


9232 


5319 


979 


979 


979 


Absorbed Flux c (ergs cm -2 s _1 ) 


5.29X10- 13 


3.62X10- 13 


5.90X10- 14 


6.50X10" 14 


8.21xl0~ 14 


Unabsorbed Flux c (ergs cm -2 s~ 1 ) 


8.66X10- 12 


4.70X10- 12 


7.52X10- 13 


1.60X10- 12 


4.34X10- 13 


Luminosity (ergs s _1 ) 


9.96 xlO 33 


5.41 xlO 33 


8.65x10 s2 


1.84X10 33 


4.99X10 32 



a All quoted error bounds are 90% confidence intervals. The best model is marked in bold where a few models are presented for the same spectra set. 
b Dcfincd as f n e rifj dV: here, d is the distance to CTB 1 (in cm) and n e and njj are the electron and H densities (in cm -3 ). 
c For the energy range 0.5-5.0 keV. 

d The numbers of degrees of freedom range from 101 to 104. 
°The numbers in brackets are the fits with fixed solar abundnaccs. 



Table 4. Summary of Fits to GIS and SIS Spectra of Southwest and Northeast Regions of CTB 1 



Parameter 


GIS2/3+SIS0/1 


GIS2/3+SIS0/1 


GIS2/3+SIS0/1 


GIS2/3+SIS0/1 


GIS2/3+SIS0/1 


GIS2/3+SIS0/1 


GIS2/3 


Region 
Model 

X 2 /DOF 

N H (x 10 22 cm- 2 ) 

kT aoft (keV) 

O b 

Nc 

Fe 

ti (10 11 cm" 3 s) 

EMi c /(47rd 2 /10- 14 ) (cm" 5 ) 


Southwest 
PHABSx 
VAPEC 

1.339 [1.33] h 
0.59 [ 0.51] 
0.30 [0.30] 
1.7 [=1] 
1.6 [1.2] 
0.5 [0.4] 

2xl0" 2 


Southwest 
PHABSx 
VNEI 

1.29 
0.68 
0.41 

1.7 

1.6 

0.1 
2 

2xl0" 2 


Southwest 
PHABSx 
(VAPEC+ 
POWER LAW) 
1.15 [1.15]9 

0.28±°f 4 
1.7 
1.6 
0. 1 

2xl0" 2 


Southwest 
PHABSx 
(VNEI+ 
POWER LAW) 

1.129 [l.l9] h 

0.57t° .tl [°- 59 ] 
0.27+q 23 [0.27] 

1.7 [=1] 
2.7(>2.0) [=1] 
0.8(>0.4) [=1] 
10 (>0.2) [15] 

2xl0" 2 


Southwest 
PHABSx 
(VAPEC+ 
VAPEC) 

1.139 
66+ 09 

U - DD -0.13 
25+ - 05 

1.7 

2.1±0.4 

0.7±0.4 

3xl0" 2 


Southwest 
PHABSx 
(VNEI+ 
VNEI) 

1.129 [i.i2] h 
0.80lp 3 [0.9] 
0.23t° g [0.22] 
1.7 [=1] 
2.0±?;5 [1.8] 

i-2±B;| [1.2] 

0.4(<5) [0.4] 
8xl0" 2 


Northeast 
PHABSx 
(VAPEC+ 
POWER LAW) 

1.19* [1.16] h 
0.47±0.32 [0.60] 
0.22+g'f 4 [0.24] 
=1.7 [=1] 
=1.6 [=1] 
=0.4 [=1] 

2xl0" 2 


r or kT hard 

t 2 (10 11 cm" 3 s) 

EM 2 C /(47rd 2 /10- 14 ) (cm" 5 ) 

Normalization 






9 o+U.4 
z - z -l.l 

4xl0" 4 


3.oli;S [3-0] 
lxlO- 3 


3-3^ 
lxlO" 3 


3.1(> 2) [3] 

o.7±5:2 [0.7] 

8xl0" 4 


1.4±0.9 

1 

2xl0- 4 g 


Absorbed Flux f 
Unabsorbed Flux f 
Luminosity 


3.52X10" 12 
3.37X10- 11 
3.88 xlO 34 


3.64xl0" 12 
5.35X10- 11 
6.16X10 34 


4.44X10" 12 
3. 34x10-" 
4.19X10 34 


4.31X10" 12 
3.50X10- 11 
4.03 xlO 34 


4.38X10" 12 
4.88X10- 11 
5.61 xlO 34 


4.43X10" 12 
l.lOxlO- 10 
1.27X10 35 


4.37X10" 12 | 
3. 02x10-" 
3.47X10 34 



a AU quoted error bounds are 90% confidence intervals and the best model is marked in bold. 

b The abundance for oxygen was frozen to 1.7 and 1 [numbers in brackets] for all fits. Since the errors are similar for both cases, only one set of errors are given. 
c Emission measure, defined as J n e nji dV: here, d is the distance to CTB 1 (in cm) and n e and nn are the electron and H densities (in cm— 3 ). 
d Thc units for the thermal component fcThard are keV. 

c Normalization for the power law component in units of photons keV -1 cm -2 s _1 at 1 keV. 

For the energy range 0.6—10.0 keV. The units for the flux are ergs cm -2 s — 1 and the units for the luminosity are ergs s — 1 . The luminosity estimates are for an assumed 
distance of 3.1 kpc. 

g,l The degree of freedom for Southwestern spectra ( 9 ) is ~1039 with both GIS and SIS data, for Northeastern spectra 1 is 394 with only GIS data. 
h The numbers in brackets are the fits with fixed solar abundnaccs. 
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Table 5. Summary of the Physical Properties of HB21 and CTB 1 



Property 


HB21 


CTB 1 


Angular size (arcmin) 


120 x 90 


34 


Distance (kpc) 


1.7 


3.1 


Physical Diameter (pc) 


59x45 


31 


Average Radius (pc) 


26 


15.3 


N H (10 21 cm" 2 ) 


0.23-0.31 


0.57-0.80 


kT or r 


0.63-0.68 keV 


0.28 keV, T = 2-3 






yuL o ive v j 


[O/Fe] 




4.3+i° 5 2 


[Ne/Fe] 




4.0j£° 


n e (cm -3 ) 


0.06 




M x h (M ) 


23.4 




X-ray pressure P/k(K cm" 3 ) 


0.9xl0 6 


2.8xl0 6 


Shock velocity (optical) (km/s) 




400 


Age (optical) (yr) 




1.6xl0 4 


Shock velocity (HI shell) d (km/s) 


124 


107? 


Age (HI shell) d (yr) 


4.5 xlO 4 


4.4xl0 4 


n (HI shell) d (cm -3 ) 


2.7 




E (HI shell) d (ergs) 


1.6xl0 51 




AM (HI shell) d (M ) 


150±10 


12±8 



a See Sections [5] and [6J 

b In the case of HB21, the listed X-ray emitting mass is for the volume 
of the whole SNR. The X-ray emitting mass of just the area of HB21 
observed by ASCA is 2.6 M . 



Fesen et all (119971 ). 



nKoo Sz Heilesl ( 199ll ). who assumed distances of 1.1 kpc and 3.2 kpc 
to HB21 and CTB 1, respectively. 
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Fig. 1. — Mosaicked broadband ASCA GIS image of HB21: the emission has been smoothed with 
a V FWHM Gaussian. The intensity range is - 7.39 xl(T 5 counts sec 1 arcmin 2 . The contours 
represent radio emission from HB21 as observed at 408 MHz with the Canadian Galactic Plane 
Survey; they range from 0.445 Jy to 0.890 Jy per 3f6 x 3f6) in steps of 0.0445 Jy/beam. The 
blue boxes and circles represent the approximate fields of view of the ASCA SIS and the GIS, 
respectively. 
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Fig. 2. — Mosaicked ROSAT PSPC image of HB21: we have superposed the same radio contours 
as in Figure [H The intensity range of this image is 10 -4 - 2.1xl0~ 2 counts sec -1 arcmin -2 . 
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Fig. 3. — Mosaicked soft band (E<1 keV - top) and hard band (E>1 keV - bottom) GIS images 
of the central region of HB21: both images have been smoothed with a V FWHM Gaussian. The 
coverages of the GIS images are marked in Figure [TJ The intensity ranges of the two images are 
- 1.02xl0" 4 and - 2.58xl0~ 5 counts sec 1 arcmin 2 , respectively. 
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Fig. 4. — GIS2 (red), GIS3 (black), SISO (green) and SIS1 (blue) spectra of the northwestern 
portion of HB21 superposed on PHABSxVNEI models (see Table [2]). 
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1 1.2 1.4 1.6 1.8 

Si abundance 

Fig. 5. — Confidence contours for the abundances of Si and S for the PHABSxVAPEC fit to the 
extracted GIS2, GIS3, SISO and SIS1 spectra of the northwestern portion of HB21. The confidence 
contours are at the la, 2a and 3a levels (see Table [2|). 



- 36 - 



o 



10:00.0 



63:00:00.0 



50:00.0 - 



40:00.0 



30:00.0 



2B:Q0.D - 



10:00.0 ■ 



62:00:00.0 



T 1 1 1 1 1 - 




M:00.0 



0:00:00.0 
HA £02000.0} 



23:56:00.0 



Fig. 6. — Mosaicked broadband ASCA GIS image of CTB 1 (color): the emission has been smoothed 
with a 1' FWHM Gaussian. The intensity range is - 5.93 xlO -5 counts sec -1 arcmin -2 . The 
contours represent radio emission as observed at 408 MHz with the Canadian Galactic Plane Survey 
(CGPS): the contours range from 0.024 Jy to 0.063 Jy per 0'9 x 0f9 in steps of 0.0035 Jy/beam. 
The blue circles depict the approximate fields of view of the ASCA GISs while the red squares 
depict the approximate fields of view of the Chandra ACIS-I chips. 
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Fig. 7. — ASCA GIS image of CTB 1 with CGPS radio contours overlaid: these radio contours are 
at the same levels as those shown in Figure [H 
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Fig. 8. — Mosaicked soft band (E<1 keV - top) and hard band (E>1 keV - bottom) GIS images 
of CTB 1: both images have been smoothed with a 1' FWHM Gaussian. The intensity ranges are 
- 6.42 xlCT 5 and - 2.85 xlCT 5 counts sec 1 arcmin 2 , respectively. The discrete source 1WGA 
J0001.4+6229 (located at RA (J2000.0) 00 h 01 m 25 s . 5, Dec (J2000.0) +62° 29' 40") becomes more 
prominent at higher energies: this source is discussed in detail in Section SJ 
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Fig. 9. — Ha image of CTB 1 in greyscale (courtesy of Robert Fesen) with blue ASCA GIS contours 
overlaid. The contour levels are 2.05, 3.04, 4.02 and 5.01 x 10~ 5 counts sec 1 arcmin 2 . Notice 
how the X-ray emission extends through the prominent gap of optical emission. 
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Fig. 10. — Three-color exposure-corrected and smoothed Chandra image of the interior of CTB 
1. This image has been made using monoenergetic exposure maps. Red, green and blue colors 
correspond to soft (0.5-1.0 keV), medium (1.0-2.0 keV) and hard (2.0-8.0 keV) emission, respectively. 
The three regions (in green) of spectral extraction are indicated: clockwise from upper left, these 
are the regions denoted as "hard," (II chip), "soft" (13 chip) and "diffuse" (12 chip) excluding 
point sources (see Section f 3 . 2 [) . The boundaries of the Chandra chips are marked with thin boxes 
(in green). Note how clumpy features with different spectral properties can be distinguished by 
noticeable differences in their colors. See Section [3. 2i 
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Fig. 11. — Chandra spectra of the three different regions of CTB 1. The black, red and green 
lines correspond to spectra from the "diffuse," "soft" and "hard" regions, respectively. The spectra 
have been fit with a PHABSxVAPEC model with the temperature of 0.28, 0.28 and 0.66 keV, 
respectively. The detailed parameters for these fits are given in Table El 



12 3 12 3 

abundance Ne abundance 

Fig. 12. — Confidence contours of O and Fe abundances and Ne and Fe abundances, respec- 
tively, showing enhanced O and Ne abudances and lower Fe abundance relative to solar, for the 
PHABSxVAPEC fit to the spectrum extracted for the "diffuse" region. The confidence contours 
are at the la, 2a and 3a levels. See Table 



0.2 0.4 0.6 0.8 1 

nH (1 22 ) 

Fig. 13. — Confidence contours for the line-of-sight column density N h and the temperature kT 
for the PHABSxVAPEC spectral fit (with fixed solar abundances) to the spectrum of the "hard" 
region (see Table [3|). Notice the bimodality for components with temperatures of kT ~ 0.28 keV 
and kT ~ 0.66 keV: this may indicate a presence of an additional thermal component besides the 
thermal component identified in the fit of the spectrum of the "diffuse" region. The confidence 
contours are at the la, 2a and 3a levels. See Section 1331 
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Fig. 14. — ASCA GIS (top) and SIS (bottom) spectra of the southwestern region of CTB 1: the 
spectra have been fit using the PHABSx(VAPEC+Power Law) model (see Table HJ). 
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Fig. 15. — ASCA GIS2 spectra of the southwest (in black with crosses) and northeast (in red with 
squares) regions of CTB 1. Both of these spectra have been fit with a PHABS x (VAPEC+POWER 
LAW) model. The parameters of the fits to these spectra are listed in Table HI 
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Fig. 16.— (left) The GIS2/GIS3 spectra of the hard source 1WGA J0001.4+6229 as fit with an 
absorbed power-law model, (right) Plot of confidence contours for Photon Index versus Njj for this 
fit: a low column density is implied, indicating that the source is Galactic rather than extragalactic. 
The confidence contours are at the la, 2a and 3a levels. See Section HJ 



